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Abstract 

We examine the growth of the Weyl curvature in two examples of 
naked singularity formation in spherical gravitational collapse - dust 
and the Vaidya spacetime. We find that the Weyl scalar diverges 
along outgoing radial null geodesies as they meet the naked singularity 
in the past. The implications of this result for the Weyl curvature 
hypothesis are discussed. We mention the possibility that although 
classical general relativity admits naked singularity solutions arising 
from gravitational collapse, the second law of thermodynamics could 
forbid their occurrence in nature. The method can also be used to 
compare the relative importance of initial data and that of the energy- 
momentum tensor in deciding the metric solution in any general case. 
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1 Introduction 



As Penrose has emphasized JT|] , a fundamental explanation for the second law 
of thermodynamics can be had only after we first understand why the entropy 
of the Universe was extremely low in the beginning. There is good cosmolog- 
ical evidence that in the very early epochs, the matter distribution itself was 
in a thermal equilibrium state with maximal entropy. Hence the matter fields 
by themselves cannot provide an explanation for the low initial entropy. The 
low entropy constraint has to come from the spacetime geometry - a suitably 
constructed geometrical quantity plays the role of gravitational entropy. This 
gravitational entropy has to be initially so very small that the net sum of the 
matter and gravitational entropy takes a value far less than it could possibly 
have. As the Universe evolves, the gravitational entropy increases as a result 
of clumping of matter and formation of irregularities. The overall entropy of 
the Universe also increases. 

A plausible geometrical quantity which could play the role of gravita- 
tional entropy is the Weyl curvature tensor [p]]. This is because at the initial 
Big Bang singularity, the Ricci curvature necessarily diverges, whereas the 
Weyl curvature is zero, if the Universe is assumed to be Friedmannian. The 
development of inhomogeneities during evolution leads to an increase in the 
Weyl curvature. Ultimately, it is expected that inhomogeneities lead to the 
formation of gravitational singularities. If cosmic censorship holds, these 
final singularities will be black-holes, and as a result of the continuous in- 
crease of the Weyl tensor, we can expect that the Weyl tensor will diverge at 
the black-hole singularity. Thus the structure of the final singularity is very 
different from the initial singularity, in so far as the behaviour of the Weyl 
curvature is concerned. 

In this way of looking at things, the validity of the second law requires 
that at the initial cosmological singularity, the Weyl curvature be zero, or at 
least ignorable in comparison with the Ricci curvature. This is the essence of 
Penrose's Weyl curvature hypothesis. Various important investigations have 
been carried out in order to formulate a precise version of the hypothesis 
[§]. For natural reasons, these investigations are largely concerned with the 
initial cosmological singularity. Here we are interested in the Weyl hypothesis 
from a different viewpoint, namely the possible relevance of this hypothesis 
for naked singularities occurring in gravitational collapse of compact objects. 
A naked singularity should be regarded as a 'T.I.F.' as well as a 'T.I. P.', 
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because both ingoing and outgoing geodesies terminate at such a singularity. 
In one particular version of the hypothesis ||, Penrose enquires if given a 
naked singularity, the Weyl curvature would diverge along ingoing geodesies, 
and go to zero along outgoing geodesies terminating at the singularity. In 
other words, one is asking if the Weyl hypothesis that is usually applied to 
an initial cosmological singularity would also hold for a naked singularity 
forming in collapse. It is this version of the hypothesis that we examine in 
our brief note. Considering the similarity between a naked singularity and 
an initial cosmological singularity in some, though certainly not all, respects 
this version of the hypothesis deserves further examination. 

This issue becomes interesting in the light of examples of naked singulari- 
ties that have been found in recent years. Here we analyse the evolution of the 
Weyl tensor in the Tolman-Bondi spacetime, and in the Vaidya spacetime, 
both of which admit naked singularities for certain initial data. As is known, 
in both cases, it is only the central shell-focusing singularity which can be 
naked. For defmiteness, we assume that it is the Weyl scalar C = C abcd C a b c d 
which has to be zero at the initial singularity. We find that the Weyl scalar 
remains zero at the center for all epochs prior to the formation of the sin- 
gularity. (The singularity is defined by the divergence of the Kretschmann 
scalar). Exactly at the (naked) singularity, however, the scalar blows up, 
both along ingoing as well as outgoing null geodesies. This appears to sug- 
gest that the version of the Weyl hypothesis under consideration does not 
hold if such naked singularities occur. 

A knowledge of the evolution of Weyl curvature is of interest also for 
reasons other than the Weyl hypothesis. Various criteria for judging the se- 
riousness of a naked singularity have been proposed - these include curvature 
strength and stability of the Cauchy horizon, among others. In addition to 
these, the divergence of the Weyl curvature at the naked singularity could 
imply that the singularity is to be taken seriously, because the Weyl tensor is 
not related to the local matter distribution, and is more directly associated 
with that part of the gravitational field which is determined by the matter 
non-locally. On the other hand, if the Weyl scalar were to be zero at the 
naked singularity, one could possibly conclude that these singularities are 
associated with the matter distribution. 

The paper is organised as follows. In Section II we study the evolution 
of the Weyl tensor in the Tolman-Bondi dust collapse. We also use this 
opportunity to present a much simplified demonstration of the occurrence 
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of the naked singularity. From this demonstration one can see that while 
the early pioneering works on inhomogeneous dust collapse established the 
formation of a naked singularity through difficult calculations, we now have 
elementary ways of rederiving those early results. In Section III the Weyl 
calculation is repeated for the Vaidya spacetime. The possible implications 
of these results are discussed in Section IV. 



2 Weyl tensor for dust collapse 

The Tolman-Bondi metric for dust collapse in comoving coordinates (t, r, 8, (f>) 
is 

ds 2 = dt 2 - Y^J^ dr 2 ~ R2 ( r > f ) dn2 W 

where R(t,r) is the area radius, and the free-function f{r) labels bound, 
marginally bound and unbound models, depending on whether f(r) is nega- 
tive, zero or positive, respectively. The field equations are 

#=/M + ^, (2) 

where dot and prime indicate partial derivative with respect to t and r, 
respectively. The quantity F(r) is equal to two times the mass inside the 
sphere labelled r. The metric solution in terms of the parameters above yields 
a singularity in the spacetime when R(r, t) = 0, where the Kretschmann 
scalar diverges. 

As has been shown earlier [|J, the central shell- focusing singularity is 
naked if the radial null-geodesic equation 

X= lim ^LtL-M (3) 
r^o.fl^o ar a - x 1 V X J 

admits one or more positive real roots X . This equation is written by 
first eliminating the variable t in favour of R, and then eliminating R in 
favour of X = R/r a , where X is the tangent to a possible outgoing geodesic. 
The constant a > 1 is chosen so that i?'/r Q_1 has a unique, finite limit as 
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R — > 0,r — > 0. The constant Ao is the limiting value of A(r) = F(r)/r a as 
r -»• 0. 

At this stage the calculations in earlier papers ||, |5[] proceed by giving a 
detailed expression for R' - this is the most involved part of the calculation. 
The simplification we present here is in the calculation of R' . Consider first 
the marginally bound case, for which the solution is 

^3/2 = r 3/2 _ l^T) t . (4) 

An initial scaling R — r at the starting epoch t — of the collapse has 
been assumed. From this equation, evaluate R' and substitute for t from 
(f|). In the resulting expression, substitute R = Xr a , and perform a Taylor 
expansion of F(r) around r = 0, so as to retain the leading non- vanishing 
term. We get that near r = 0, 

U ' X + J =r ?+3/2-3«/2 (5) 



r a-l ^ x 

Here, 9 = —qF q /3F , and q is defined such that in a series expansion of the 
initial density p(r) near the center, the first non-vanishing derivative is the 
gth one (=p q ), and F q = p q / (q+S)q\. Now a unique definition of a is obtained 
by setting the power of r as zero in the second term, giving a = 1 + 2q/3. 
This then reproduces the result of the R' calculation performed earlier || , in 
a simpler manner. From this point on, the naked singularity analysis through 
evaluation of roots of @ proceeds as before, but the overall calculation is 
now more economic. 

For the non-marginally bound case, the solution of the Tolman-Bondi 
equations is 



R 3/2 G(-fR/F) = r 3/2 G(-fr/F) - yJF{r)t (6) 

where G is a known elementary function ||. R' can be evaluated as before, 
and then we eliminate t and substitute R = Xr a . The power-series expansion 
for f(r) near r = is of the form 

f(r) = f 2 r 2 + f 3 r 3 + Ur 4 + ... (7) 

which implies that the argument (—f'R/F) of G on the left-side of (^) goes 
to zero as r — >• 0. The derivative of G(—p) w.r.t. its argument p = rf/F is 
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obtained by differentiating (El), which gives 

dG{-p) _ 3G _ 1 
d(— p) 2p py/1 + p' 

Using this, one can now perform a Taylor-expansion of G, F and / about 
r = to get exactly the same expression for R' as given in @ above, where 
9 is now given by 



e = q i 



G(-h/F ) 



2 i 
Fo 



3/g+2 
2/2 



1 + 



A" 

2F 



+ 



g+2 
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The constant g is now defined such that the first non- vanishing derivative of 
the initial density is the gth one, and/or the first non- vanishing term in the 
expansion for f(r), beyond the quadratic term, is of order r q+2 . The constant 
a is again equal to (1 + 2g/3). Thus the R' calculation is simplified for the 
non-marginal case as well. 

We return now to the calculation of the Weyl tensor, which is defined as 

Cabcd = Rabcd + 9a[dRc\b + 9b[cRd]a + -Rga[c9d]b, 

where square brackets denote antisymmetrization. The nonvanishing inde- 
pendent components for the Tolman-Bondi metric are 



The Weyl scalar is given by 



_F_ 

1+/ \R 3 
F | F' 
2R ~r 6R' 



F'R' 



ri R'* l f 

W010 — -r—F I — 

C2020 

C*2121 

C3030 

C*3131 

C*3232 



+ 



3R 2 R' 



F(r)R n 



6(1+/) 1 2R(l+f) 

C2020 sin 2 9 
C2121 sin 2 ^ 



(-FR+^) sin* 9 



(9) 



At the start of evolution at t = 
see that the Weyl scalar at r 




(10) 



0, where the scaling R = r has been chosen, we 
= is zero (the mass-function F(r) grows as r 3 



6 



near r = 0.) Using the above Tolman-Bondi solution we find that the scalar 
remains zero at r = throughout the non-singular phase of the evolution. 
However, at the epoch of naked singularity formation, we know that in the 
neighborhood of the singularity, we have along an outgoing geodesic the 
relation R = X Q r a , X being the finite tangent. Also, R' is given by fl5|), 
with X = Xq. With these substitutions, we see that in the limit of approach 
to the singularity, the scalar is equal to 

C(X 0> r) = - 12F ° 2g ° 2 -,r^~ a \ (11) 

XI (X + flo/v^) 

For those initial data that lead to a naked singularity, the range of a is 
5/3 < a < 3. Hence we see by letting r go to zero that the Weyl scalar di- 
verges at the naked singularity along outgoing geodesies. In so far as ingoing 
geodesies are concerned, these exist for all initial data, i.e. a takes all values 
greater than unity. Hence the Weyl scalar necessarily diverges along ingoing 
geodesies too. 

The effect of the energy-momentum tensor on the geometry can be found 
by evaluating the scalar R a pR a/3 and comparing it to the Weyl scalar. We 
get in the limit of approach to the naked singularity, 

TP' 2 Q P 2 

B „B aP - — - — 2 r 6{1 ~ a) (12) 

R R X$(x + 6 /VX^) 

We find that this scalar diverges at the same rate as the Weyl scalar, and 
the latter does not dominate over the former. Notice however that unlike 
the Weyl scalar, the invariant R a pR a/3 is non-zero at the initial epoch, and 
evolves with time. Finally we note that for the points r ^ the Weyl scalar 
is in general non-zero initially - it then evolves with time and diverges at the 
singularity R(t, r) = 0. We recall that the singularity at r ^ is necessarily 
covered. 
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3 Weyl tensor for Vaidya spacetime 



The collapse of null dust is described by the Vaidya metric 

ds2 = ^ _ 2n^)\ ^ 2 _ _ r , dn2 (i3) 

where v is the advanced time, and m(v) is the mass-function. The energy- 
momentum tensor is 

1 dm 

T ik = euiu k = - — ^—UiU k (14) 

where = —5^ - this represents the radial inflow of radiation into an ini- 
tially flat region of spacetime. The spacetime outside the radiation region is 
Schwarzschild. In the self-similar model, which is the case that we consider 
here, the mass-function in the cloud is given by 2m(v) = Xv, i.e. it is linear. 
A curvature singularity forms when a shell hits the origin r = 0. It is known 
that the singularity at r = 0, v = is naked for A < 1/8 |J. 

The non-zero independent components of the Weyl tensor are given by 

j — y _ 2m(v) 



1010 — 

2m(v) 
r 



r 1 — m ( v ) ( i 

(-^2020 — — U — 



/^i _ m(v) 

<--2120 — — ~ 

C3232 = —2m(v)r sin 2 9 
CVnso = C2020 sin 2 9 



(15) 



^3030 — <^2020 

C3130 = C2120 sin 9 
The Weyl scalar is given by 



. 48m 2 12AV 
C(v,r) = — ^ = — j-. (16) 



Prior to the formation of the singularity, the scalar is zero at the inner edge 
v = 0, r 7^ of the cloud. Naturally, it is also zero in the flat region to 
the interior of the collapsing cloud. When the naked singularity forms, the 
finite tangent along an outgoing geodesic from the singularity is defined as 
X = v/r (see M). Thus we get 



12A 2 X 2 



C(X ,r) = ^-f°- (17) 
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where Xo is the limiting value of the tangent. The Weyl scalar diverges at 
the naked singularity. The scalar R a pR al3 is however zero, and in this case, 
the Weyl scalar dominates. 

4 Discussion 

Our aims in this paper were twofold. Firstly, to compute the Weyl tensor 
in the approach to the naked singularity in the Tolman-Bondi and Vaidya 
spacetimes, and secondly to consider what the results could mean for the 
Weyl curvature hypothesis. We find that in these models, the Weyl scalar 
behaves in a peculiar manner, remaining zero at the origin throughout the 
non-singular phase of the evolution, and then suddenly jumping to infinity 
along outgoing geodesies, at the epoch of singularity formation. The same is 
true along ingoing geodesies. In so far as the Weyl hypothesis is concerned, 
it appears reasonable that the hypothesis apply in the same manner to cos- 
mological initial singularities and to naked singularities ||. In the examples 
considered here, we find that the curvature structure of the T.I.F.s is sim- 
ilar to that of the T.I.P.s, which does not support the hypothesis. Also, 
since naked singularities are expected to arise in strongly inhomogeneous 
situations, for which the Weyl tensor is typically high, one could speculate 
that the Weyl scalar should generically diverge in the past along outgoing 
geodesies, rather than go to zero. 

If we do assume that the Weyl hypothesis must be valid, then these results 
offer an interesting possibility, namely that the second law of thermodynamics 
could forbid the occurrence of naked singularities in nature. If the Weyl 
hypothesis applies both to cosmological initial singularities and to naked 
singularities arising from gravitational collapse, then validity of the second 
law requires that only those naked singularities can actually occur in nature 
for which the Weyl scalar goes to zero along outgoing geodesies. The second 
law would hence disallow those naked singular solutions which disagree with 
the Weyl hypothesis. If the Weyl scalar diverges generically along outgoing 
geodesies, the second law would generically forbid naked singularities from 
occurring in nature. 

Thus while classical general relativity admits naked singularities in grav- 
itational collapse, their occurrence in the real world could be thermodynam- 
ically forbidden. The situation will then be analogous to the advanced wave 
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solutions in classical electrodynamics - solutions that are permitted by the 
theory but thermodynamically disallowed in the real world. The cosmic cen- 
sorship hypothesis would then hold, though not in classical general relativity 
as such, but in the real universe, through the intervention of thermodynam- 
ics. However, this would mean attributing a fundamental status to the Weyl 
Curvature Hypothesis, which could possibly involve quantum gravity. 

Although it was the notion of connection of the Weyl tensor with non- 
uniformity of energy momentum distribution which led Penrose to think of 
its relation to entropy, it must be kept in mind that the Weyl tensor has 
a certain part which is determined by the 'initial' conditions of the field 
equations. These are the metric and its derivatives specified on a chosen 
space-like surface. The Weyl scalar going to zero from TIF would mean that 
one would have to have unique initial conditions on a partial Cauchy surface 
chosen near the singularity formation event, as it approaches the ideal point. 
Thus, Penrose seems to indicate that there was only a unique choice of the 
metric solution || , presumably at a limiting space-like surface near the ideal 
point. 

Also, given the fact that the divergence of the Weyl scalar and R a pR al3 
are comparable, if not the Weyl contribution dominating, one would be led 
to think that whenever a naked singularity occurs, the role played by the 
initial conditions in deciding the metric of the universe is as important as 
that of the energy momentum tensor, if not more. 

It must be noted that the conjecture includes all singularities generically 
and not just the cosmological ones. We find that it is not true for Tolman 
Bondi dust and Vaidya radiation collapse scenarios. 

This leaves the issue of the factors deciding the initial conditions open. 
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